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INTRODUCTION 
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Advances  in  fields  of  scientific  endeavor  are  typically  made  in  response  to  new 
performance  requirements.  This  certainly  has  been  the  case  in  the  area  of  computational 
plasticity  where  there  have  been  increasing  demands  from  government  and  Industrial  or- 
ganizations for  analytically  determining  accurate  stress,  strain,  and  displacement 
fields  in  a wide  range  of  complex  structures  subjected  to  severe  environmental  emd 
mechanical  loading  conditions.  These  considerations  have  been  a motivating  force  be- 
hind the  development  of  general  purpose  programs  for  the  plastic  analysis  of  structures. 
The  "Decade  of  Progress"  in  computational  plasticity  is  represented  by  the  general 
availability  and  acceptance  of  these  multi-purposb  computer  programs  for  nonlinear 
structural  analysis.  This  achievement  has  been  the  direct  result  of  advances  in  the 
areas  of  structural  mechanics  and  computer  sciences.  Specifically,  the  burgeoning  de- 
velopments in  finite  element  methods,  made  possible  by  advances  in  computing  hardware 
and  software,  formed  the  foundation  from  which  efficient  and  accurate  algorithms  for 
plastic  analysis  could  be  developed. 

Several  alternative  algorithms  have  evolved  and  are  currently  being  used  in  pro- 
grams for  plastic  analysis.  Each  of  these  has  computational  advantages  and  liabilities 
so  that  the  analyst  is  presented  with  choices  among  acceptable  procedures. 

These  procedures  are  reviewed  to  facilitate  a rational  selection  among  competing 
procedures  that  best  suit  the  needs  of  a particular  problem  class. 
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The  development  of  «.l((orlthms  vlthin  the  freaevorh  of  finite  element  computational 
techniques  for  plastic  analysis  was  accomplished  independently  from  the  development  of 
appropriate  constitutive  relations.  Initially,  only  the  simplest  plasticity  theories 
were  used.  With  the  availability  of  a computational  tool  that  could  treat  problems  of 
a more  complex  nature  than  previously  practical  (within  constraints  of  computer  time) 
analysts  could  turn  their  attention  towards  the  development  of  plasticity  theories  that 
could  more  accurately  predict  the  essential  features  of  experimentally  observed  behavior 
without  the  penalties  of  introducing  unrealistic  simplifications. 

Implicit  in  the  development  of  any  plasticity  model  are  assumptions  associated  with 
the  behavior  of  the  material . The  number  and  degree  of  these  assusiptlons  effect  the  gen- 
erallty  of  the  resulting  model  and  its  compatibility  with  actual  material  behavior.  Some 
assumptions  generally  employed  are  listed  with  a discussion  of  their  implications.  This 
Is  followed  by  a synoptic  presentation  of  several  models  that  incorporate  these  assumptions 
or  attempt  to  overcome  their  restrictions. 

The  models  discussed  have  either  been  incorporated  into  available  general  purpose 
plastic  analysis  programs  or  have  been  proposed  for  possible  Implementation. 

ASSVJMi'TlONS  OF  PUiSTlC  BEHAVIOR 

1.  The  existence  of  an  initial  yield  condition  defining  the  elastic  limit  of  the 
material  in  a multlaxlal  stress  state.  This  assumption  is  most  often  used  and  represents 
a mathematical  convenience  that  facilitates  computational  efficiency.  Its  applicability 
varies  with  the  material  under  consideration. 

The  twvi  popular  yield  criteria  for  structural  materials  are  those  attributed  to  von 
Mises  and  Tresca.  The  former  implies  that  yielding  begins  at  an  arbitrary  material  point 
whenever  the  combination  of  stresses  is  such  that  the  strain  energy  of  distortion  per 
imit  volume  at  the  point  is  equal  to  the  correspondittg  energy  developed  in  a bar  unlax- 
ially  stressed  to  the  elastic  limit. 
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The  Treaca  condition  states  that  inelastic  action  at  any  point  in  a body  begins  only 
when  the  maxifflum  shearing  stress  on  soae  plane  through  the  point  reaches  a value  equal  to 
the  maximuin  sheai'ing  stress  in  a tension  specimen. 

2.  The  existence  of  loading  surfaces  that  define  the  limits  of  elastic  and  plastic 

behavior  beyond  initial  yielding  the  response  after  initial  yielding  differs  among 

various  plasticity  theories.  This  post  yielding  response,  called  the  hardening  rule,  is 
described  by  specifying  a subsequent  yield  surface,  termed  the  "loading  surface,"  which 
is  a convenient  mathematlceQ  idealisation  of  some  macroscopically  observed  behavior.  The 
"consistency"  condition  requires  that  the  stress  state  at  any  point  remain  on  the  loading 
surface. 

3.  Plastic  strain  rates  are  llneai'ly  related  to  their  corresponding  stress  rates  by 

means  of  a flow  rule  (the  basis  of  constitutive  relations  in  the  treatment  of  plasticity). 
A flow  rule  that  is  generally  used  to  describe  elastlc-ldeally  plastic  behavior  is  the 
Prandtl-Reuss  relation,  which  is  a generalization  of  the  Levy-Mlses  equations.  The 
Pi'andtl-Reuss  assumption  is  that  the  plastic  strain  Increments,  proportional  to 

the  corresponding  stress  deviator,  where  the  instantaneous  non-negative  value  of  the 

constant  of  proportionality  is  left  to  the  Inventiveness  of  the  user  of  these  equations. 
The  concept  of  an  effective  stress  or  effective  plastic  strain  (generatlly  in  terms  of  one 
or  a combination  of  their  corresponding  Inviurlants)  is  in  Itself  an  assumption  that  is 
usually  Introduced  to  reduce  the  complexity  of  a multiaocial  situation  to  one  that  cam  be 
related  to  uniaxial  behavior  [l].  Thus,  the  proportionality  parameter  can  be  the  ratio  of 
the  effective  plastic  strain  Increment  to  the  prevailing  effective  stress. 

A more  general  approach  to  determining  a flow  rule  is  the  use  of  the  concept  of  a 
plastic  potential.  The  assumption  is  made  that  there  exists  a scalar  function  of  stress, 
say  from  which  the  component  of  plastic  strain  Increments  may  be  determined  (e.g., 

proportional  to  3f/30jj).  If  f(Oj^j)  represents  the  yield  surface  in  stress  space,  then 
the  above  assumption  represents  a result  of  Drucker’s  postulate  (2],  which  states  that 
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th»  vork  done  by  an  external  agency  during  a complete  cycle  of  loading  and  unloading  must 
be  non-negative.  Furthermore,  this  assumption  leads  to  an  incremental,  or  associated 
linear  flow  theory  of  plasticity,  in  which  the  increment  of  plastic  strain  (strain  rate) 
is  in  the  direction  of  the  outward  normal  to  the  surface  represented  by  stress 

space,  at  the  current  value  of  stress.  A strain  rate  vector  deviating  from  the  outward 
normal  to  the  yield  surface  In  a direction  Independent  of  the  stress  rate  vector  consti- 
tutes non-assoclated  flow  theories  of  plasticity.  Non-associated  flow  theories  are  par- 
ticularly suitable  for  work-softening  materials  and  can  reasonably  fit  the  behavior  ob- 
served in  some  soils.  Associated  and  non-associated  flow  theories  are,  in  general,  dis- 
tinct from  the  deformation  theory  of  plasticity,  in  which  the  total  plastic  strains  are 
related  to  the  final  stress  state.  According  to  this  latter  theory,  a relationship  between 

final  states  of  stress  and  strain  exists  for  any  given  loading  process  unloading  being 

specified  by  a separate  law. 

1*.  Materials  are  Isotropic  with  respect  to  initial  yielding  — whether  Induced  by 
previous  cold-working  or  as  a result  of  an  anisotropic  distribution  of  crystals,  most 
structural  materials  exhibit  some  form  of  Initial  anisotropy.  The  von  Mlses  and  Tresca 
criteria  assume  isotropy  with  respect  to  the  orientation  of  the  stresses  and  their  sense 
(tensile  or  compressive). 

5.  Plastic  Incompressibility  Incompressibility  assumptions  are  generally  em- 

ployed in  moat  nonlinear  structural  analysis  programs.  Investigations  of  the  Influence 
of  hydrostatic  pressure  on  the  plastic  response  of  metals  have  been  considered  by  Brld^nnan 
(3).  Despite  some  evidence  to  the  contrary  [4],  most  investigations  suggest  that  hydro- 
static pressure  has  little  or  no  effect  on  the  initial  and  subsequent  behavior  of  metals 
under  quasi-static  loading  rates.  This  has  been  most  recently  confirmed  by  Fung,  et  al. 


The  use  of  plastic  Incompressibility  requires  treating  a variable  Poisson  ratio, 
l.e.  an  elastic  value,  and  a vaJ.ue  equal  to  one-half  In  the  plastic  range.  The  impli- 


cations of  this  are  discussed  In  assumption  6. 

The  treatment  of  the  plastic  yielding  of  nonmetals,  such  as  clay,  ice,  and  concrete, 
represents  an  area,  of  Investigation  in  which  the  effects  of  hydrostatic  pressure  are 
significant.  The  plasticity  theories  employed  for  the  treatment  of  these  materlad.s  are 
non-assoclated  linear  flow  types. 

6.  Yielding  and  subsequent  response  are  Insensitive  to  rate  of  deformation  — the 
nature  of  the  constitutive  equations  of  plasticity  for  rate-sensitive  material  behavior 
has  been  the  subject  of  several  notable  investigations  (see  survey  paper  of  [6]).  These 
investigations  have  shown  that  limits  on  dislocation  velocity  and  rate  of  dislocation  den- 
sity for  most  structural  materials  are  responsible  for  the  generally  observed  fact  that 
plastic  flow  Is  retarded  with  Increasing  strain  rate. 

7.  Elastic  unloading  and  coincidence  of  yield  and  loading  surfaces  — this  is  eua 
idealization  of  material  behavior,  that  implies  that  \inloading  from  some  plastic  state  to 
neighboring  state  results  in  a change  in  the  elastic  state  only.  Furthermore,  implicit 
in  most  analyses  eire  the  assumptions  that  additional  plastic  strains  can  occur  only  upon 
reloading  to  a stress  state  beyond  that  of  Initial  unloeullng  and  subsequent  behavior 
(assuming  no  reverse  yielding  has  occurred)  is  identicaLl  to  that  which  would  have  been 
obtained  had  unloading  never  occurred. 

8.  The  total  kinematic  strain  may  be  dec(nposed  into  elastic  and  plastic  components 

this  assumption  is  enployed  in  the  linear  flow  theories  of  plasticity  and  is  used  in 

developing  many  of  the  incremental  constitutive  equations.  The  assimptlon  is  mathematic- 
ally convenient  when  used  in  conjunction  with  the  concept  of  a well-defined  yield  surface 
and  is  general ly  valid  for  small  deformations.  In  addition,  this  assumption  facilitates 
the  use  of  a veurlable  Poisson  ratio  when  plastic  incompressibility  is  enployed. 
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In  addition  to  the  above  assumptions,  most  successful  applications  of  classical  small 
strain  flow  theories  of  plasticity  have  been  generally  limited  to  situations  in  which  the 
loading  is  monotonically  Increasing  and  the  ratio  of  the  various  load  components  are  held 
constant  (proportional  loading).  Applications  of  flow  theory  to  problems  involving  severe 
changes  in  stress  ratio  amotxg  the  vai'lous  components  of  stress,  such  as  occurs  during 
plastic  buckling,  have  been  shown  to  be  unconservative  because  the  reduction  in  the  in-plane 
shear  stiffjiess  is  generally  underestimated.  Hutchinson  has  discussed  this  problem  exten- 
sively in  (7l.  Deformation  theory,  on  the  other  hand,  has  been  shown  to  provide  consis- 
tently better  correlation  with  results  for  bifurcation  problems  for  plates  and  shells, 

[8,  9).  This  situation  Illustrates  a significant  fact  that,  with  few  exceptions,  appears 
to  have  been  overlooked  in  many  experimental  investigations.  While  the  sise,  shape,  and 
definition  of  yield  and  subsequent  loading  surfaces  are  of  considerable  importance,  equal 
emphasis  should  be  given  to  a proper  description  of  the  plastic  hardening  modulus  under  a 
variety  of  conditions. 

I'yclic  loading  situations  involving  reversed  plastic  flow  during  which  the  material 
response  may  exhibit  cyclic  hardening  or  softening  behavior  falls  into  the  category  of 
problems  of  general  interest  in  which  classical  plasticity  theories  have  met  with  limited 
success.  Applications  of  cyclic  loading  conditions  to  several  problems  ranging  from 
membrane  stressed  sheets  to  plates  and  shells  have  been  presented  in  [lO,  ll].  In  these 
studies,  "stabilized"  material  behavior  was  assumed. 

Finally,  most  applications  of  plastic  analysis  have  been  limited  to  isothermal  condi- 
tions. The  influence  of  elevated  temperatures  on  material  response  (elastic  and  Inelastic) 
is  generally  treated  within  the  framework  of  theories  of  viscous  strain  (creep).  For  a 
great  many  problem  areas,  ranging  from  metal  forming  processes  to  the  analysis  of  nuclear 
reactor  components,  elevated  stress  eind  temperature  levels  of  short  duration  must  be  tol- 
erated. In  these  problems  plastic  strain  development  may  be  quite  significant.  Thus, 


non-isotheraal  plasticity  theory  is  required.  Furthermore,  it  is  desirable  that  such  a 
theory  include  the  effects  of  temperature  on  the  ^-lastic , as  well  as  the  plastic  proper- 


ties of  the  materials. 


MODELS  OF  PLASTIC  BEHAVIOR 


The  following  is  a brief  description  of  some  models  for  plastic  behavior.  Several 
of  these  models  have  been  developed  within  the  framework  of  the  assumptions  mentioned; 
others  have  been  developed  specifically  to  overcome  one  or  more  of  the  restrictions  as- 
sociated with  these  assumptions.  The  primary  function  of  the  models  is  to  provide  a set 
of  constitutive  relations  that  cem  be  used  to  describe  the  response  history  of  nonlinear ly 
deforming  media. 

ISOTROPIC  HARDENING 

This  theory,  proposed  by  Hill  [12]  and  Hodge  [13],  assumes  that  during  plastic  flow 
the  loading  surface  expands  uniformly  about  the  origin  in  stress  space,  maintaining  the 
same  shape,  center,  and  orientation  as  the  yield  surface.  Figure  1 illustrates,  on  the 
basis  of  a simplification  to  a two-dimensional  plot,  the  yield  and  loading  surfaces  when 
the  stress  state  shifts  from  point  1 to  2.  Unloading  and  subsequent  reloading  in  the  re- 
verse direction  will  result  in  yielding  at  the  stress  state  represented  by  point  3.  The 
path  2-3  will  be  elastic,  and  0-2  is  equal  to  0-3. 

Isotropic  representation  of  work  hardening  does  not  account  for  the  Bauschinger  effect 
exhibited  by  most  structural  materials.  In  fact,  contrary  to  observations,  this  theory 
implies  that,  because  of  work  hardening,  the  material  will  exhibit  em  increase  in  the 
compressive  yield  stress  equal  to  the  increase  in  the  tensile  yield  stress.  Furthermore, 
since  plastic  deformation  is  an  anisotropic  process,  it  cannot  be  expected  that  a theory 
that  predicts  isotropy  in  the  plastic  range  will  lead  to  realistic  results  when  complex 
loading  paths,  involving  changes  in  direction  of  the  stress  vector  in  stress  space  (not 
necessarily  completely  reversed),  are  considered. 


SLIP  THEORY 


Utilizing  the  physical  concept  of  slip  surfac  » in  crystals,  Batdorf  and  Budiansky 
[ll*]  have  developed  a theory  that  describes  a loadinj?  surface  that  is  distorted  relative 
to  the  yield  surface  and  previous  loading  surfaces.  This  theory  predicts  the  fonaation 
of  corners  at  the  instantaneous  stress  state  on  the  loading  surface  during  plastic  deform- 
ation. After  some  prestrain  the  yield  locus  is  the  minimum  surface  through  the  point  of 
prestrain  and  the  initial  yield  locus.  A representation  of  the  growth  of  the  yield  func- 
tion in  going  from  a stress  state  at  the  origin  0 to  the  final  state  represented  by  point 
3 is  given  in  Fig.  2.  In  this  figure,  the  unshaded  region  is  that  enclosed  by  the  yield 
surface,  and  the  various  shaded  regions  indicate  the  stages  (I,  II  and  III)  in  the  forma- 
tion of  the  loading  surfaces  in  going  from  stress  state  0 to  3. 

Since  the  stress  state  is  almost  always  in  a corner,  the  resulting  constitutive  rela- 
tion between  stresses  and  strains  becomes  quite  complex.  For  this  reason  this  theory  is 
rarely  selected  for  application. 

PIECEWISE  LINEAR  PLASTICITY 

In  this  representation,  the  yield  surface  consists  of  a finite  number  of  plane  surfaces 
whose  intersections  constitute  comers.  The  oldest  and  most  widely  used  piecewise  linear 
yield  surface  is  that  associated  with  the  Tresca  yield  condition.  The  loading  surface  is 
assumed  also  to  consist  of  pleu^e  surfaces,  and  the  subsequent  hardening  behavior  can  be 
classified  as: 

1.  The  hardening  rule  of  independent  plane  loading  surface  one  of  the  earliest 

discussions  of  this  representation  of  the  hardening  behavior  is  given  in  [15]  and  is  illus- 
trated in  Fig.  3(A).  As  seen  from  this  figure,  in  which  and  are  the  only  nonzero 
stress  components,  a loading  path,  0-1-2,  in  any  quadrant  of  the  stress  plane  does  not 
affect  the  loading  surface  in  the  remaining  quadrants.  Thus,  this  hardening  rule  does 
not  take  the  Bauschinger  effect  into  account. 


2.  The  hardening  rule  of  interdependent  loading  surface  this  type  of  hardening 

rule,  origineilly  proposed  by  Hodge  [l5]t  is  a generedization  of  the  hardening  rxxle  des- 
cribed in  the  previous  paragraph.  By  specifying  a dependence  between  the  plames  that 
compose  the  loading  surface,  a loading  path  intersecting  any  one  plane  of  this  surface 
may  effect  changes  in  each  of  the  remaining  planes.  As  illustrated  in  Fig.  3(B),  this 
hardening  rule  can  be  used  to  specify  any  piecewise  linear  loading  surface  and  is  capable 
of  taking  the  Bauschinger  effect  into  account. 

A special  case  of  the  Interdependent  loading  surfaces  is  considered  in  [l6].  It  is 
assumed  that  plastic  strain  is  due  to  slipping  along  three  Independent  slip  planes,  along 
any  one  of  which  the  shear  is  a maximum.  Piecewise  linear  stress-strain  relations  are 
written  in  terms  of  coefficients  representing  the  hardening  behavior  of  the  material. 

These  coefficients  are  functions  of  stress  etnd  ewe  dependent  upon  a lineeur  strain-harden- 
ing rule  employed  in  the  analysis.  By  specifying  the  correspondence  between  various  seg- 
ments of  the  yield  surface  and  the  slip  planes,  toteLL  plastic  strains  for  eujy  loading  ewe 
computed  as  the  sum  of  the  contributions  from  the  three  independent  sets  of  slip  plames. 

It  is  further  assumed  that  the  corresponding  segments  of  the  yield  surface  must  maintain 
a constant  elastic  reuige  from  positive  to  negative  yielding.  An  illustration  of  the  sub- 
sequent loading  surfaces  determined  in  this  way  is  shown  in  Fig.  3(C).  It  is  seen  from 
this  figure  that  the  Bauschinger  effect  can  be  tedcen  into  a'.  :ount. 

A comprehensive  review  of  the  piecewise  linear  strain-hardening  theory  of  plasticity 
is  presented  in  [l?]. 

KINEMATIC  HARDENING 

The  hardening  behavior  postulated  in  this  theory  assumes  that  during  plastic  deforma- 
tion the  loading  surface  translates  as  a rigid  body  in  stress  space,  maintaining  the  size, 
shape,  and  orientation  of  the  yield  surface.  The  primary  aim  of  this  theory,  due  to 
Prager  [l8,  I9],  is  to  provide  a means  of  accounting  for  the  Bauschinger  effect. 
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For  piecewise  linear  yield  surfaces,  kinematic  heurdenlng  may  be  considered  to  be  a 
special  case  of  the  hardening  rule  of  Interdependent  loading  surfaces.  However,  it  Is 
not  limited  to  piecewise  linear  yield  surfaces. 

An  illustration  of  klnonatlc  heurdenlng,  as  applied  in  conjunction  with  the  von  Nises 
yield  curve  in  the  <^,  plane,  is  provided  in  Fig,  1*.  The  yield  surface  and  loading 

surface  are  shown  in  this  figure  for  a shift  of  the  stress  state  from  point  1 to  point  2. 
The  translation  of  the  center  of  the  yield  sxirface  is  denoted  by  a... 

As  a consequence  of  assuming  a rigid  translation  of  the  loading  surface,  kinematic 
hardening  predicts  an  ideal  Bauschinger  effect  for  completely  reversed  loading  conditions. 
A modification  to  this  theory,  proposed  by  Ziegler  [20],  eliminates  Inconsistencies  that 
arise  when  Praiger's  original  model  Is  used  In  a subspace  of  stress. 

Although  originally  devised  to  be  used  in  conjunction  with  llnecur  strain  heurdening 
behavior,  this  model  has  been  used  for  materials  exhibiting  nonlinear  hardening  [lO,  ll], 
and  has  been  further  generalized  to  cyclic  loading  involving  work-hardening  and  work- 
softening behavior  [21]. 

A model  of  combined  kinematic  and  isotropic  hardening  in  which  the  subsequent  loading 
surfaces  expand  and  translate  is  presented  by  Hodge  [22]. 

MECHANICAL  SUBLAYER  AND  FIELDS  OF  WORKHARDENING  MODULI  MODELS 

A technique  to  model  the  arbitrary  nonlinear  mechanlcsil  behavior  of  a solid  by  means 
of  a parallel  assemblage  of  elastic  ideally  plastic  solids  can  be  traced  to  Duwez  [23], 
with  extensions  by  White  [2h]  and  Bessellng  [25].  This  modeling  concept  equates  the  in- 
tegrated effect  of  a network  of  ideedly  plastic  solids  to  the  actual  behavior.  An  exten- 
sion of  this  model,  proposed  by  Mroz  [26,  27]  to  account  for  the  work-hardening  behavior 
of  metals  under  cyclic  loading  conditions.  Introduces  the  notion  of  a field  of  work-hard- 
ening moduli  and  the  variation  of  this  field  during  the  course  of  plastic  deformation. 

I 
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In  this  proposed  model,  a stress-strain  curve  of  an  Initially  isotropic  material  is 
represented  by  n linear  segments  of  constant  tangent  plastic  moduli,  as  shown  in  Fig.  5. 


j In  stress  space,  this  approximation  can  be  represented  by  n hypersurfaces  f^,  f^^ 

where  f^  is  the  initial  yield  surfaces,  and  f^^  to  f^  define  regions  of  constant  work- 
hardening  modul i . 

. Figure  6 illustrates  these  hypersurfaces  in  the  plane  for  an  initially  iso- 

* tropic  material.  As  seen  in  this  figure,  the  surfaces  f^,  f^  ...  f^  are  similar  and  con- 

j centric,  and  for  simplicity  are  schematically  represented  by  a family  of  circles.  If  we 

consider  porportlonal  loading  in  the  direction,  corresponding  to  o in  Fig.  6,  and  if 
I we  assume  that  the  surfaces  can  experience  a rigid  translation  without  experiencing  a 

change  of  size  or  orientation,  then  when  the  stress  state  reaches  point  A on  Fig.  6,  the 
I surface  f^  will  translate  until  it  reaches  the  circle  f^  at  the  stress  corresponding  to 
j point  B.  Tiie  circles  f^  and  f^^  translate  together  until  point  C is  reached,  where  now 

f , f and  f^  are  attached  at  a common  point  of  contact.  For  unloading  and  subsequent 

I reversed  loading,  when  the  stress  reaches  a point  corresponding  to  point  E (Fig.  6),  re- 
verse plastic  flow  occurs  and  the  surface  f^  translates  downward  along  the  until 

I it  reaches  the  8urf..ce  f^  at  F.  Mroz  further  proposes  that  the  curve  of  reverse  loading 
in  Fig.  B .loin  the  curve  OA’B'O  that  is  obtained  by  symmetry  with  respect  to  the  origin 
j from  OABC.  Thus,  the  curve  of  reverse  loading  EFG  is  uniquely  defined  by  the  curve  of 
primary  loading,  represented  by  an  equation  of  the  form  o » f(e).  If  a new  coordinate 
I system  (o,  e),  with  origin  at  C is  used,  we  have  for  the  curve  CEFO,  0/2  ■ f(e/?). 

I ’ This  relation  is  usually  referred  to  as  the  Masing  relation  [28],  and  is  a useful  rule 

for  describing  steady  cyclic  behavior. 

In  the  generalization  of  this  model  to  nonproportional  loading,  it  is  assvuned  that 
during  translation  of  the  hypersurfaces  the  individual  surfaces  do  not  inters*  ^ but  con- 
I secutively  contact  and  push  each  other.  It  should  be  noted  that  when  fj^  tends  to  infinity. 


in  which  case  the  work-hardening  modulus  is  constant  (the  work-hcurdenlng  curve  being 
represented  by  a straight  line),  the  theory  proposed  by  Mroz  is  identical  to  Prager's 
klnsaiatic  hardening  model. 

The  further  generalization  of  the  theory  of  work-hardening  moduli  is  associated  with 
an  expansion  or  contraction  of  the  surfaces  f^,  f^^  f^  so  that  transitory  phenomena 

(work-stiffening,  work-softening,  or  non-isothermal  conditions)  can  be  treated.  Thus, 
the  hypersurfaces  f,.  are  not  constants  but  functions  of  a monotonlced.ly  increasing  scalar 
parameter  during  plastic  flow.  One  suggestion  for  the  scalar  is  presented  in  (27). 

TWO  SURFACE  THEORIES 

Eisenberg  and  Phillips  [29]  presented  an  early  generalization  of  conventional  plas- 
ticity theories  to  account  for  the  phenomenon  of  noncolncident  yield  and  reloading  stress 
states.  To  account  for  this  behavior  a two-surface  plasticity  theory  was  proposed  involv- 
ing a yield  surface  completely  enclosed  by  a prescribed  locuiing  surface.  The  site  of  the 
yield  surface  always  remains  unchanged.  The  loading  surface  on  the  other  hand  varies  in 
size,  coinciding  with  the  yield  surface  for  elastic  behavior  and  containing  the  stress 
point  for  loading  beyond  initial  yielding.  During  unloading  the  yield  surface  remains 
unchanged  and  the  loading  surface  shrinks  to  acconmodate  the  stress  state  until  it  coin- 
cides with  the  yield  surface.  Upon  unloading  cmd  subsequent  reloading  the  two  surfaces 
separate. 

The  Incorporation  of  this  behavior,  should  it  be  a significant  factor  for  a material 
under  consideration,  does  not  appear  to  pose  any  difficulties  beyond  those  normally  asso- 
ciated with  conventional  theories. 

A comprehensive  and  satisfying  generalization  of  the  concept  of  a two-surface  plas- 
ticity theory  was  proposed  by  Dafallas  and  Popov  [30].  In  this  theory  the  concept  of  a 
bounding  surface  is  introduced.  This  surface,  always  enclosing  the  yield  and  subsequent 
loading  surface  in  stress  space,  is  used  to  model  complex  loading  situations,  including 


cyclic  loading  involving  hardening  and  softening  behavior.  This  work  appears  to  have 


been  motivated  by  the  general  observation  of  the  noncoincidence  of  the  yield  and  loading 


stresses  previously  discussed.  The  approach  to  multiaxlal  loading  as  presented  by  Dafcdias 


and  Popov  is  not  tied  to  any  hardening  law  and  appears  to  be  sufficiently  general  and 


sound  to  warrant  further  examination 


A modification  of  Mroz's  fields  of  work-hardening  moduli  model  has  been  recently 


proposed  by  Krieg  [3l].  This  modification,  similar  to  the  model  proposed  by  Dafalias 


and  Popov  [30],  replaces  all  but  two  of  the  discrete  surfaces  specified  by  a Mroz  model 


by  a continuum  of  Intermediate  loading  surfaces  whose  distribution  is  prescribed.  The 


two  surfaces,  are  represented  by  an  inner  curve  labelled  by  Krieg  as  the  loading  surface. 


and  an  outer  curve,  termed  the  limit  siirface.  These  two  surfaces  separate  the  material 


behavior  into  three  distinct  zones:  an  elastic  zone  contained  within  the  loading  surface 


an  asymptotic  plastic  zone  outside  the  limit  surface,  and  a so-celled  "metaelastic"  zone  be- 


tween the  two  surfaces.  On  the  basis  of  a uniaxial  stress-strain  curve,  these  zones  are 


Joined  by  a continuous  function,  generated  for  a variety  of  situations  including  reversed 


loading.  Both  the  loading  and  limit  surface  can  vary  according  to  a combined  kinematic 


isotropic  heu*denlng  behavior.  The  motion  of  the  loading  surface  is  Identical  to  that 


assumed  by  Mroz.  For  the  general  multiaxlal  case  the  theory  requires  the  retention  of 


three  vectors  and  three  scalars,  a small  increase  over  the  two  vectors  required  for  kine- 


matic hardening  alone 


ANISOTROPIC  THEORIES 


One  of  the  first  treatments  of  the  plastic  flow  of  an  initially  anisotropic  metal 


was  suggested  by  Hill  [12].  In  this  theory  an  orthotropic  yield  criterion  was  assumed  to 


be  quadratic  in  the  stress  components,  and  to  reduce  to  the  von  Mises  law  when  the  degree 


of  anisotropy  was  small.  In  Cartesian  coordinates  the  function,  suggested  by  Hill,  tcUies 


the  following  form 


'■<'>»>  ■ "ijVj  - " • ° 

j 

where  the  contracted  tensor  notation  is  used,  i,  J , k > 1,  2,  . . . 6,  and  the  parameter 
i N^j  are  constants  related  to  the  six  yield  stresses  in  the  principal  directions  of  aniso- 

tropy. A consequence  of  this  assumed  function  is  an  initial  rotation  of  the  von  Mises 
ellipse  in  the  space. 

Simileu:  anisotropic  theories  have  been  suggested  by  severed  investigators,  most 
notably  Jackson,  Smith  and  Lankford  [32],  Dom  [33]  and  Hu  [3**].  Implicit  in  Hill's 

i theory  ore  the  following  assumptions: 

i 

I Orthotropic  anisotropy 

1 The  principal  axes  of  anisotropy  either  coincide  with  the  principal,  stress  axes  and 

t the  principal  strain  axes  or  the  transformation  between  the  aoces  are  known 

The  principal  axes  of  anisotropy  do  not  rotate  during  plastic  flow 


No  distinction  is  made  between  tensile  and  compressive  stresses 

The  anisotropic  coefficients  remain  unchanged  during  plastic  flow. 

A generalization  of  Hill's  equation  for  anisotropic  plasticity,  and  one  that  combines 
Isotropic  emd  kinematic  hardening  is  presented  by  Boltov  and  Sawczuk  [33]>  In  compact 
notation  this  theory  assumes  a yield  condition  in  the  following  form 

‘■'“k’  ■ "ij  ‘“i  - “j)  - 1 ■ 0 (2) 

Unlike  previous  investigations,  the  fourth  order  tensor  of  anisotropic  coefficient,  , 
is  a prescribed  function  of  the  plastic  strains,  and  hence  changes  during  the  course  of 
plastic  flow.  The  kinematic  hardening  peurameter,  a^,  is  also  a function  of  the  plastic 
strains.  In  a comparison  of  experimental  results  for  combined  tension  and  torsion  it 


was  found  that  the  proposed  theory  better  suited  the  experimental  data  than  the  conven- 
tional kinematic  hardening  theory. 


More  recently,  several  investigations  have  been  concerned  with  defining  a strength 
criterion  for  orthotropic  materials  with  specific  reference  to  advanced  filamentary  com- 
posites. Generally,  a criterion  of  this  type  is  used  to  define  a "failure  surface"  in 
stress  space.  These  failure  surfaces  may  be  useful  for  defining  yield  surfaces  of  ortho- 
tropic  metals.  Notable  among  these  proposed  theories  is  the  work  of  Tsai  and  Wu  [36]. 

The  basic  assumption  associated  with  this  anisotropic  strength  criterion  is  a failure 
surface  in  the  following  form 


-1.0 


The  parameters  and  N^j  are  strength  tensors  of  the  second  and  fourth  rank  respectively 
and  are  functions  of  an  appropriate  number  of  Independent  material  strengths.  The  linear 
term  in  the  above  equation  can  be  used  to  account  for  the  difference  between  tensile  and 
compressive  induced  "failure". 

RATE-SENSITIVE  MODELS  ^ 

The  survey  paper  of  Lee  [6]  lists  some  200  references  associated  with  the  investiga- 
tion and  application  of  dynamic  plasticity.  Various  functional  representations  and  varying 
degrees  of  complexity  have  been  proposed.  Many  of  these  are  impractical  for  use,  either 
because  of  their  extreme  complexity  or  because  they  have  been  simplified  beyond  the  point 
of  usefulness. 

The  most  popular  constitutive  relation,  suggested  by  Malvern  (37l.  employs  the  con- 
cept of  a reference,  or  static,  stress-strain  function.  The  dynamic  stress-strain  behavior 
if  determined  from  the  static  curve  in  some  prescribed  manner  as  a function  of  the  strain 
rate.  For  uniaxial  conditions  this  relation  is  in  the  following  form 


- 1 ♦ 
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where  n and  D are  material  constants,  0^  and  o are  the  static  and  dynaaic  yield  stresses, 
respectively,  and  4 is  the  strain  rate. 

A generalization  of  the  above  actuation,  as  suggested  by  Pertyna  [38],  replaces  the 
uniaxial  strain  rate  component  by  the  tensor  invariant  of  the  plastic  deformation  rate 

Concepts  associated  with  tine- independent  plastic  models  may  be  extended  within  the 
framework  of  the  above  assumptions  to  treat  rate-sensitive  yielding.  On  this  basis  the 
rate-dependent  loading  function  may  be  written  in  the  following  form 


‘ij' 


(6) 


where  is  the  second  stress  invariant,  and  K is  the  yield  stress  in  shear. 

Another  aspect  of  rate  sensitivity,  in  addition  to  its  effect  on  yield  stress  levels, 
is  the  variation  of  the  strain-hardening  behavior  of  the  material  with  varying  levels  of 
strain  rate.  In  conventional  flow  theories  this  dependence  must  be  Incorporated  in  the 
plastic  modulus  which  will  now  be  a function  of  strain  rate  as  well  as  stress  level. 
NON-ISOTHERMAL  MDDELS 

A treatment  of  non-lsothermal  response  of  structures  experiencing  simultaneous  changes 
in  load  and  temperature  clearly  requires  a substantial  modification  of  existing  isothermal 
procedures.  Without  consideration  for  irreversible  thermodynamics  and  for  the  uncoupled, 
quasi-static  problem,  the  treatment  of  thermoplastic Ity  requires  the  elastic-plastic 
constitutive  equations  to  account  for  the  influence  of  temperature  on  the  elastic  co- 
efficients (primarily  restricted  to  Young’s  modulus),  yield  stress,  plastic  hardening 
coefficient,  and  rate  of  thermal  expansion. 
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A theory  of  non-isothermnl  deformation  of  rigid,  work  hardening  solids  is  presented 
by  Prager  [39l-  In  this  model  a general  constitutive  law  for  plastic  flow  is  developed 
that  is  h08»geneous  of  order  one  in  the  rates  of  temperature,  stress,  and  strain.  A yield 
criterion  is  chosen  to  be  a function  of  the  state  variables  of  stress,  plastic  strain, 

p 

degree  of  hardening  and  temperatures,  l.e,  f ■ f(o,  e , h,  T).  The  motivation  for  in- 
cluding the  degree  of  hardening  (or  hardening  modulus)  term,  h,  is  to  provide  the  materiel 
with  a '’memory”  of  some  previous  plastic  deformation  that  will  affect  subsequent  behavior. 
Conditions  for  loading,  neutral  loading  and  unloading  from  some  plastic  state  are  repre- 
sented by  the  following  conditions 


15  “ * I?  O’’ 


> 0 loading 
■ 0 neutral  loading 


< 0 unloading 


If  the  problem  is  dlscretised  with  respect  to  space  and  some  monotonlcally  increasing 
parameter  that  can  be  used  to  prescribe  the  applied  mechanical  uid  thermal  history  of  load- 
ing, then  an  incremental  relation  between  plastic  strain,  stress  and  temperature  can  be 
represented  by  the  following  linear  relation 


Ae*’  - AOj  + Cat 


where  the  contracted  tensor  notation  is  used  and  C is  the  conventional  plasticity  coiiq)liance 
matrix,  modified  to  account  for  the  effect  of  temperature  on  the  hardening  characteristics 
of  material;  and  C is  a function  of  stress  and  accounts  for  the  influence  of  temperature 
varying  yield  stresses. 

Examples  of  the  treatment  of  thermoplastic  behavior  formulated  within  the  framework 
of  finite  element  methods  is  presented  by  Ueda  and  Yamakawa 
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[Ul].  Problems  in  both  geometric  and  material  nonlinearities  are  considered.  The  treat- 
ment of  changing  elastic-plastic  properties  with  temperature  is  of  particular  Interest. 
Accounting  for  varying  material  properties  is  a formidable  computational  problem  since  it 
may  lead  to  nonsymmetrlc  stiffness  influence  coefficients.  This  situation  is  discussed 
within  the  context  of  two  general  solution  procedures  presented  in  the  next  section,  e.g. 
the  initial  strain  approach  and  the  tangent  modulus  method. 

MODIXS  WITHOUT  A YIELD  SURFACE 

A theory  that  attempts  to  circumvent  the  ambiguities  associated  with  defining  a spe- 
cific yield  surface  and  prescribing  hardening  as  the  subsequent  translation,  expansion,  or 
distortion  of  that  surface  has  been  proposed  by  Valanls  [1*2,  l*3l.  In  his  theory  of  plas- 
ticity eui  intrinsic  time  parameter,  independent  of  external  clock  time,  is  chosen  to  be  a 
monotonically  increasing  function  of  deformation.  One  obtains  the  stress  response  by 
monitoring  this  history  of  strain  associated  with  this  deformation. 

From  a physical  viewpoint  the  theory  has  many  advantages.  Phenomena  such  as  cross- 
hardening, noncoincident  yield  and  loading  points  and  cyclic  hardening  are  capable  of  being 
described.  On  the  other  hand,  the  absence  of  a yield  surface  does  not  facilitate  the  com- 
putational effort  in  a general  purpose  discrete  model  analysis.  In  fact,  it  may  introduce 
additional  effort. 

ideally  plastic  behavior 

The  plasticity  models  discussed  thus  far  are  used  to  describe  the  hardening  behavior 
of  materials  subsequent  to  initial  yielding.  They  attempt  to  describe  the  process  of  haurd- 
ening  graphically  as  an  expansion,  translation,  distortion,  etc.,  of  the  initial  yield  sur- 
face. For  the  case  of  elastic-ideadly  plastic  behavior  the  yield  surface  is  assumed  to  re- 
main unchanged.  In  [ll],  the  treatment  of  multiaxial  ideally  plastic  behavior  requires 


that  the  stress  increment  vector  be  tangent  to  the  yield  surface  and  the  plastic  strain 
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increment  vector  be  normal  to  the  loading  surface.  The  conditions  on  the  stress  rate  is 
the  consistency  condition  ensuring  the  stress  state  to  remain  on  the  yield  surface.  This 
condition  provides  a linear  relationship  astong  the  various  components  of  stress  rate.  The 
condition  on  the  strain  rate  provides  s linear  relation  sBK>ng  the  various  components  of 
plastic  strain  rates.  The  combination  of  independent  increments  of  stress  and  plastic 
strain  are  subsequently  determined  as  a function  of  the  increment  of  total  strain  as  out- 
lined in  [111. 

REVIEW  OF  STATIC  SOLUTION  PROCEDURES 

In  the  following  it  is  assumed  that  the  reader  is  sufficiently  familiar  with  the 
derivation  of  the  necessary  components  of  elastic  stiffness  influence  coefficients,  form- 
ulated within  the  framework  of  the  displacement  method  of  finite  element  analysis.  The 
matrix  equation  relating  generalised  displacements  and  loads  may  be  written  in  the  follow- 
ing matrix  notation 

Ku  - P (9) 

where  K is  the  matrix  of  conventiorsd.  elastic  stiffness  influence  coefficients,  u is  the 
vector  of  generalised  displacmnents  euid  P is  the  vector  of  genered.ised  loads. 

. Basically,  the  procediire  used  to  solve  the  small  displacement  plasticity  problem  may 

• be  divided  into  two  categories:  in  one  the  effects  of  plutlclty  are  accounted  for  directly 

in  the  stiffness  Influence  coefficients;  the  second  treats  plasticity  as  an  effective  load 
that  is  used  in  conjunction  with  the  applied  mechanical  and  thermal  loads  for  general  equil- 
ibrium. The  latter  is  referred  to  as  the  residual  force  or  Initial  strain  method,  and  the 
former  la  termed  the  tangent  modulus  method.  A derivation  of  their  corresponding  equations 
from  virtual  work  principals,  the  relationship  between  the  iq>proaches,  and  a discussion  of 
several  variations  are  presented  in  [lO],  [ll]  and  [1»1*]. 
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The  governing  equation  associated  with  each  technique  ia  customarily  written  in  incre- 
mental form.  The  tangent  modulus  procedure,  for  example,  uses  the  following  form 

K.p  Au  - AP  (10) 

where  the  prefix  notation  A denotes  the  rate  of  change  of  the  parameter  following  it  with 
respect  to  any  monotonlcally  increasing  function  of  time,  and  represents  the  matrix  of 
elastic-plastic  stiffness  Influence  coefficients.  The  plasticity  model  chosen  for  use  is 
explicitly  contained  in  K.J,. 

The  residual  force  method  uses  the  following  matrix  relation 

KAu  - AP  ♦ AQ  (11) 

where  the  vector  Ai^  is  the  plastic  residual  load  vector.  For  the  initial  strain  approach 
Aij  is  the  product  of  the  initial  strain  matrix  and  the  Incronents  of  plastic  strain;  AQ  ■ 
K^ACp.  For  the  initial  stress  approach,  AQ  is  the  product  of  the  initial  stress  matrix  and 
the  Increments  of  plastic  strain.  The  plasticity  model  explicitly  enters  the  analysis 
through  the  calculation  of  the  increment  in  plastic  strain. 

It  should  be  mentioned  that  the  treatment  of  Ideally  plastic  behavior  may  be  incorpor- 
ated in  the  tangent  modulus  method  or  both  approaches  associated  with  the  residual  force 
method. 

The  most  obvious  way  of  solving  either  set  of  equations  for  a complete  history  of 
response  is  a forward  Euler  integration  scheme  in  which  the  applied  load  history  is  divided 
into  incremental  load  steps.  Assuming  the  solution  is  known  at  acme  arbitrary  step  1,  we 
desire  the  solution  for  the  increment  of  displacement  associated  with  the  applied  increment 
of  load.  The  displacement  solution  is  used  to  obtain  total  strain  increments  (from  kinematic 
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considerations).  Using  the  elastic  and  plastic  constitutive  equations  and  the  assumption 
that  allows  the  decomposition  of  the  total  strain  into  elastic  and  plastic  components,  we 
can  determine  Increments  of  plastic  strain  and  stress.  This  information  is  used  as  a basis 
for  going  on  to  the  next  step,  and  the  process  is  repeated  until  the  desired  history  has 
been  obtained.  Thus,  the  solution  to  the  nonlinear  response  is  obtained  as  a sequence  of 
linear  solutions  in  which  the  set  of  stiffness  influence  coefficients  are  modified  by  a 
residual  force  to  maintain  equilibrium. 

% 

The  advantages  associated  with  each  procedure  tire  contained  in  the  following  synopsis; 

o The  tangent  modulus  method  permits  the  use  of  larger  load  Increments  than  the 
residual  force  method  because  no  approximation  need  be  made  for  the  current 
increment  of  plastic  strain.  Step  site  is  thus  limited  by  the  approximations 
inherent  in  the  assumptions  of  flow  theory  of  plasticity. 

o The  residual  force  method  retains  the  same  set  of  elastic  stiffness  coefficients 
throughout  the  analysis  so  that  calculations  associated  with  forming  the  stiff- 
ness matrix  and  solving  the  governing  equations  need  be  performed  only  once. 
Furthermore,  this  method  may  be  used  for  the  analysis  of  materials  that  exhibit 
strain-softening  behavior  for  materials  that  require  a non-assoclative  flow 
rule  to  develop  their  plastic  constitutive  relations,  and  for  non-lsothermal 
plasticity  problems  where  the  material  compliances  vary  with  temperature.  In 

each  of  these  latter  situations  use  of  the  tangent  modulus  approach  results  in 

* 

a nonsymmetric  stiffness  matrix  that  is  generally  costlier  to  deal  with  than  the 
symmetric  case. 

The  disadvantages  associated  with  each  procedure  are  the  mirror  images  of  the  above 
advantages;  namely: 

o The  tangent  modulus  method  requires  a successive  reformulation  and  redecom- 
position of  the  stiffness  coefficient  matrix,  K^.  For  optimum  computational 
efficiency  this  set  of  coefficients  should  be  positive  definite,  thus  pre- 
cluding materials  that  exhibit  strain  softening,  or  have  non-assoclative  flow 

rules,  or  consider  plastic  material  properties  that  vary  with  temperature. 
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To  maintain  stability  and  accuracy,  the  residual  force  method  requires  small 
load  steps  to  be  used  in  the  successive  linearisation  procedure,  primarily 
because  a predicted  value  of  the  plastic  strain  increments  is  used  in  Eq.  (ll). 


The  size  of  the  load  increment  is  an  obvious  queintity  that  is  commonly  used  aS  a 
basis  of  comparison  between  competing  successive  linearization  schemes.  However,  it  is 
not  and  should  not  be  construed  as  the  overriding  factor.  For  example,  it  is  generally 
recognized  that  in  using  the  tangent  modulus  approach  greater  load  increments  than  those 
required  for  a comparable  solution  using  the  residual  force  method  are  possible.  On  the 

I 

surface,  this  would  represent  a strong  argument  favoring  the  former.  It  is  the  experience 
of  this  writer,  however,  that  the  ratio  of  computational  times  required  for  a full  decom- 
position and  subsequent  solution  (as  required  in  the  tangent  modulus  approach)  to  the 
time  required  to  obtain  the  solution  alone  (as  in  the  residual,  force  method)  can  be  great 
enough  to  offset  the  gains  associated  with  the  larger  load  increment.  The  question  of 
which  approach  is  more  cost-effective  is  now  not  quite  so  clear. 

With  either  approach  we  can  expect  the  successive  linearization  procedure  to  drift  frcan 
a true  equilibriiun  position  for  the  nonlinear  response.  This  drifting  is  a combined  result 
of  truncation,  the  successive  linearization  procedure  and  the  fact  that  information  not  yet 
available  is  required  for  a true  solution.  In  the  case  of  the  tangent  modulus  method,  this 


information  is  necessary  to  form  the  updated  stiffness  coefficients  K^;  for  the  residues, 
force  method,  the  value  of  the  plastic  strain  increment  is  required  to  form  AQ.  Several 
techniques  (iterative  and  non-iterative)  can  be  applied  to  reduce  the  amount  of  drifting  to 
tolerable  levels.  Many  of  these  techniques  are  reviewed  by  Tillerson  [U5].  With  respect 
to  both  the  tangent  modulus  and  residual  force  method  the  simplest  corrective  procedure 
involves  the  introduction  of  an  equilibrium  correction  term  that  may  be  added  as  a load 
vector  at  regular  intervals  (not  necessaurily  each  step)  in  the  incrementaO.  procedure.  ) 

The  equilibrium  correction  term  at  the  1th  step  is  represented  by  in  the  following  ■ 

for  the  tangent  modulus  approach  ' 

- f^  (12)  » 

. i 
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where  f represents  the  internal  forces  at  the  start  of  the  step.  For  the  residual  force 


method 


= Ku^  _ 


which  can  be  shown  to  be 


= K»  (Ae^  - Ae^"^) 
P P 


If  the  system  is  in  equilibrium  at  the  ith  step  R is  a null  vector.  Otherwise,  it  must 
be  added  to  the  incremental  equation  at  the  i + 1 step.  This  procedure  represents  a one- 
step  iteration.  It  can  eLLso  be  used  in  an  iterative  manner  within  a load  step  until  the 
correction  term  converges  to  a predetermined  small  value. 

Several  other  techniques  have  been  suggested  and  successfully  incorporated  to  improve 
the  accuracy  associated  with  an  incremental  Euler  integration  scheme  for  the  static  plastic 
analysis  of  structures.  These  procedures  are  summarized  in  the  following: 

1.  The  use  of  mid-increment  correction  procediires:  This  technique  used  by  Felippa 
[U6]  in  conjimction  with  the  tangent  modulus  method  to  solve  problems  of  combined  geometric 
and  material  nonlinearity  requires  the  solution  of  two  sets  of  equations  within  a single 
load  step.  Half  the  incremental  load  is  applied  and  the  solution  used  to  evaluate  the 
tangent  modiilus  coefficients  at  that  mid-increment.  The  problem  is  then  resolved  using  the 
entire  increment  of  loewi  and  the  mid-increment  tangent  modulus.  A two-stage  Iterative  pro- 
cess that  allows  larger  load  Increments  to  be  used  in  conjunction  with  the  residual  force 
method  is  described  by  Vos  [**7].  This  procedure  involves  employing  an  iterative  solution 
to  the  problem  at  the  beginning  of  the  increment,  assuming  a fixed  normal  to  the  yield  sxir- 
face,  and  then  re-evaluating  the  normal  at  the  midpoint  of  the  plastic  increment. 

2.  Stress  scaling  to  ensure  consistency:  It  is  possible  that  the  successive  linear- 
ization procedure  results  in  a situation  where  at  the  end  of  an  increment  of  load  the  stress 
state  of  an  element  may  lie  outside  the  yield  surface.  Stricklin,  et  al.  [U8]  have  shown 


that  stress  scaling  back  to  the  yield  surface  to  ensure  consistency,  and  subsequently  dis- 
tributing the  unbalanced  nodal  forces  during  the  next  load  Increment  Is  a conqputatlonally 
efficient  procedure. 

3.  The  successive  implementation  of  the  constitutive  relations  on  a sublncremental 
basis  within  each  load  step  [U9]:  This  allows  for  a larger  load  step  to  be  used  to  eval- 
uate displacements  and  total  strain.  These  are  then  divided  into  a number  of  smaller  in- 
crements that  eu:e  used  with  the  constitutive  relations.  Subincrementation  is  of  particular 
Importance  when  using  the  tangent  modulus  method  where  it  is  imperative  to  limit  the  number 
of  complete  solutions  to  the  governing  equations  without  violating  the  plasticity  theory. 

h.  The  use  of  multipoint  stress  locations:  The  use  of  multiple  points  within  an  ele- 
ment for  stress  (or  strain)  evaluation  assumes  an  Important  role  in  an  elastic-plastic  anal- 
ysis. Since  the  constitutive  relations  are  imposed  at  each  of  these  points,  the  plastic 
strain  variation  and  hence  the  elastic-plastic  boundary  within  an  element  is  dependent  on 
the  number  and  location  of  the  stress  points,  the  user  in  effect  is  modeling  for  the  anti- 
cipated plastic  resi)onse. 

This  is  a particularly  useftil  notion  when  using  the  residual  force  method  since  the 
finite  element  model  can  be  formulated  based  on  the  necessary  accuracy  required  to  obtain 
the  displacement  field  and  the  definitions  of  the  number  euid  location  of  the  stress  points 
determines  the  effective  load  vector.  This  procedure  was  followed  for  a three-dimensional 
isoparametric  hexahedra  element  by  Levy  et  al.  [50]  for  several  elastic  and  elastic-plastic 
problems.  The  cost-effectiveness  of  the  procedure  is  not  as  clear  when  applied  to  the 
tangent  modulus  method  since  the  nimber  and  location  of  the  stress  points  are  an  integral 
part  of  the  evaluation  of  the  stiffness  matrix. 

5.  The  use  of  substructuring  techniques:  The  objective  of  substructxuring  procedures, 
or  static  condensation,  is  to  reduce  computational  costs  and  Improve  the  accuracy  associated 
with  a nonlinear  analysis  by  eliminating  a substantial  amount  of  computations.  The  method 
is  particularly  well  suited  to  the  initial  strain  method  for  treating  problems  of  contained 


1 

I 

i 


plastic  flow,  and  where  the  regions  of  pure  elastic  behavior  and  elastic-plastic  behavior  ! 

can  be  distlr^guished  and  estimated  with  a high  level  of  confidence.  In  this  case,  since 

the  stiffness  matrix  is  not  updated,  a reduced  matrix  obtained  by  static  condensation  is  |; 

I 

initially  formed  and  then  used  throughout  the  plastic  analysis.  Another  associated  saving  I 

l! 

is  obtained  by  recovering  stresses  and  strains  only  in  elements  in  the  reduced  set.  The 
substructuring  technique  can  also  be  used  with  additional  effort  imd  complexity  with  the 
tangent  stiffness  method  [5l].  A study  of  substructuring  as  applied  to  three  particular 
problems  Involving  small  contained  regions  of  plastic  flow  in  which  the  residual  force 
method  is  used  is  presented  in  152).  The  results  of  this  study  may  be  summarized  by  stat- 
ing that  a reduction  of  an  order  of  magnitude  in  computing  time  requirements  may  be  real- 
ized per  incremental  step  by  Incorporating  a substructuring  option. 

Finally,  the  accuracy  and  efficiency  of  any  solution  to  a nonlinear  structural  analysis 
problem  depends  upon  a number  of  other  factors,  not  the  least  of  which  is  the  user's  exper- 
ience in  setting  the  finite  elmnent  idealization  with  respect  to  size,  arrangement , and 
type  of  elonents  used.  These  factors  have  been  discussed  by  Waltz  et  al.  [53l  and  Trucke 
[5I4]  for  linear  elastic  systems  and  are  equally  applicable,  if  not  more  critical  for  elastic- 
plastic  analyses. 

NONLINEAR  DYNAMICS 

A prerequisite  to  developing  a successful  nonlinear  dynamic  analysis  capability  is  a 
sound  knowledge  of  the  options  and  pitfalls  associated  with  nonlinear  static  analysis. 

Many  of  the  comments  concerning  efficiency  and  accuracy  discussed  in  connection  with  static 
analyses  are  directly  applicable  to  nonlinear  dynamics.  The  additional  ingredients  are  the 
inclusion  of  an  Inertia  term  to  the  governing  equations  and  the  introduction  of  a time  in- 
tegration scheme  to  solve  the  equations  of  motion.  Much  attention  [55-63]  has  been  given 
to  the  development  and  evaluation  of  discrete  methods  to  numerically  integrate  the  equations 
of  motion.  Little  will  be  added  here  to  the  technical  content  of  these  references,  rather 
we  will  outline  the  various  options  and  pitfalls  that  exist  with  respect  to  nonlinear  dynamic 
analysis  and  make  some  summary  and  experience  based  comments. 


At  the  outset  It  should  he  stated  that  a significant  factor  effecting  time  step  size 
for  a dynamic  plastic  analysis  is  the  plasticity  theory.  That  is,  the  time  step  must  he 
such  that  the  assumptions  intrinsic  to  the  plasticity  theory  are  not  violated.  This  is  in 
contrast  to  linear  dynamic  antilysls  for  which  step  size  is  controlled  entirely  hy  numerical 
stability  and  accuracy  considerations. 

Basically  the  procedures  used  to  solve  the  equations  of  motion  of  hoth  lineeu*  and  non- 
linear dynamic  analysis  are  divided  into  two  categories,  direct  time  Integration  and  modal 
superposition.  Some  of  the  distinctions  of  the  two  methods  are  discussed  helow. 

DIRECT  INTEGRATION  PROCEDURES 

The  coupled  equations  of  motion  for  an  undamped  discrete  system  ceui  he  written  as 

MAii  + Af  » AP  + R (15) 


where  M is  the  mass  matrix,  Af,  AP  the  incremental  internal  and  externed  force  vectors,  re- 
spectively, and  R the  residual  load  vector.  Indices  are  omitted  from  Eq.  (15)  but  It  Is 
implied  that  the  equation  is  written  for  the  i + 1 increment,  l.e.,  in  passing  from  t^^  to 
t^  ^ The  two  unknowns  in  Eq.  (15)  are  the  incremental  vector  of  internal  forces  Af,  and 
the  incremental  acceleration  vector.  All.  The  mass  matrix  may  he  formulated  on  the  basis  of 
a "consistent"  or  lumped  mass  approach.  The  distinction  between  the  two  being  whether  the 
mass  is  represented  by  means  of  finite  element  interpolating  functions  or  is  directly  dis- 
cretized by  lumping  components  at  nodes.  The  former  approach  leads  to  a banded  matrix 
while  the  latter  leads  to  a dlagonsQ  matrix.  Which  approach  to  use  depends  on  the  integra- 
tion scheme  employed  in  the  analysis  [61(].  Two  approaches  to  the  solution  of  Eq.  (15)  can 
be  taken  each  effecting  the  form  of  the  equation  to  be  solved  in  a sequence  of  time  steps. 

In  the  first  type,  referred  to  as  explicit,  Af  is  obtained  entirely  from  previous  informa- 
tion so  that  elements  of  the  mass  matrix  are  the  only  coefficients  of  the  unknowns  (acceler- 
ations or  displacements).  In  the  second  type,  termed  implicit,  combinations  of  the  mass 
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and  stiffness  matrices  are  combined  to  form  coefficients  of  the  unknowns.  The  choice  of 


which  method  to  use  is  clearer  for  linear  problems  than  for  nonlinear  ones,  with  implicit 
methods  overwhelmingly  used  for  structural  dynamics  and  explicit  methods  for  problems, 
where  high  frequency  response  is  significant,  as  in  the  treatment  of  wave  propagation 
effects. 

Explicit  integrators  are  generally  conditionally  stable  [55]  with  the  critical  time 
step  inversely  proportional  to  the  highest  frequency  in  the  discrete  model.  Implicit  inte- 
grators are  generally  unconditionally  stable  and  tend  to  filter  out  the  higher  frequency 
response.  This  allows  for  larger  time  steps,  the  choice  of  which  is  controlled  by  the 
modes  necessary  to  predict  the  essential  features  of  the  response.  For  linear  problems 
using  a constant  time  step,  both  methods  lead  to  coefficient  matrices  that  are  constat 
through  the  entire  response  spectrum. 

The  complicating  factor  for  nonlinear  problems  is  a consequence  of  the  change  in 
stiffness  due  to  plasticity.  The  operational  choices  on  which  method  to  use  in  this  case 
are  not  unlike  the  choices  between  using  the  temgent  modulus  or  initial  strain  approaches 
for  static  nonlinear  problems  since  they  involve  trade-offs  between  smaller  less  "costly" 
time  steps  for  explicit  integration  versus  larger  but  relatively  more  "costly"  time  steps 
for  implicit  integration.  Reference  to  the  term  "costly"  here  is  related  to  the  degree  of 
complexity  and  magnitude  of  subsidiary  computations  during  each  time  step.  Some  of  the 
distinctions  of  the  two  methods  are  outlined  below. 

Explicit  Integrator 

The  explicit  formulation  solves  Eq.  (15)  directljr  as 

MAii  = AP  - Af  + R (l6) 

with  Af  obtained  from  past  information.  For  example,  if  Af  ■ K.pAu,  then  using  the  central 

2 .. 

difference  operator  Au  * 2Au^  - Au^^  At  AUj^.  Quantities  without  indices  imply  the 


1 + increment.  There  is  small  distinction  between  the  temgent  stiffness  and  initial 
strain  approaches  for  explicit  integration  since  Af  ■ K^Au  for  the  tangent  stiffness  method 
€md  Af  » KAu  + AQ  for  the  initial  strain  approach.  In  either  case  the  end  result  is  a re- 
sidual Internal  load  vector.  Further,  the  implication  is  that  the  operations  are  performed 
on  the  element  level  rather  than  performing  operations  on  the  larger  globed  syston.  Since, 
as  stated  above  the  equations  of  motion  require  only  the  incremental  vector  of  internal 
forces,  the  notion  of  a stiffness  matrix  con  be  discarded  edtogether  (as  suggested  by 
Belytchko  [65]),  with  Af  obtained  directly  from  tui  integral  Involving  stress  and  the  matrix 
that  maps  displacement  to  strains.  In  either  case  calculations  Involve  successive  solutions 
to  Eq.  (l6)  with  incremental  alterations  to  the  right  hand  side.  This  becomes  a simple 
procedure  with  a diagoneul  mass  matrix  since  the  operation  involves  simple  division.  With 
a consistent  mass  matrix  an  initial  factorization  can  be  performed  so  that  subsequent  solu- 
tions require  only  a forward  and  back  substitution  of  banded  trlanguleir  matrices.  In  pas- 
sing, we  remark  that  Kreig  and  Key  [61*]  indicate  that  there  is  an  improvement  in  accuracy 
when  the  lumped  mass  approach  is  used  with  an  explicit  Integrator  because  the  errors  in  the 
discrete  system  caused  by  this  combination  tend  to  be  coiuiterbalonclng.  Because  of  the 
ease  in  obtaining  solutions  to  Eq.  (l6),  the  computation  time  for  an  explicit  method  becomes 
strongly  dependent  on  the  element  level  stress/strain  recovery  and  the  formation  of  Af. 
Computer  costs  are  therefore  directly  tied  to  the  number  of  elements  in  the  discrete  model 
and  the  number  of  time  steps  necessary  in  the  analysis.  This  leads  to  the  major  drawback 
of  explicit  methods,  namely  that  more  refined  models  have  an  increased  frequency  spectrm, 
which  for  numerical  stability  require  a smaller  time  step.  Consequently,  there  is  a com- 
plementary effect  caused  by  a larger  set  of  elements  in  combination  with  smaller  time 
steps.  Because  of  this  situation  a break-even  point  occurs  when  the  economics  of  simpler 
calculations  are  overridden  by  the  requirement  of  ever  smaller  time  steps. 

The  most  popular  explicit  technique  currently  in  use  is  the  constant  step  central  dif- 
ference operator  [55].  It  has  been  our  preference,  however,  to  use  a variable  time  step 
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Modified  Adams  Predictor-Corrector  method  [66]  that  is  explicit  in  the  predictor  and  implicit 

I in  the  corrector  solutions.  The  advantage  of  this  method  is  that  the  time  step  is  automa- 
tically chosen  to  reflect  current  system  stiffness  and  dynamic  response.  Chu'  experience 
has  been  that  this  method  automatically  choses  time  steps  near  those  required  by  the  central 
difference  method. 

1 

Use  of  an  explicit  Integrator  fits  naturally  into  a nonlinear  analysis  since  there  are 
very  little  additional  calculations  required  compared  to  a linear  analysis.  The  smaller 
time  steps  are  consistent  with  the  initial  premise  that  the  assumptions  of  plasticity 
theory  are  not  to  be  violated.  However,  as  the  frequency  spectrum  increases,  stability 
requirements  cause  the  critical  time  step  to  be  reduced  so  that  ultimately  the  method  be- 
comes uneconomical.  We  conclude  this  discussion  with  a direction  for  future  research  that 
leads  to  an  explicit  integrator  capable  of  filtering  out  higher  frequency  response  while 

1 - 

maintaining  the  desired  accuracy  in  the  lower  frequency  regime. 

! 

I Implicit  Integrator 

The  implicit  formulation  is  based  on  difference  operators  that  contain  both  the  cur- 

j 

i . rently  unknown  acceleration  and  displacement.  When  an  operator  of  this  form  is  substituted 

, into  Eq.  (15)  this  yields 


KAu  • AP  + R (17) 

M 

where  R “ K.p  + gZt^’  ^d  ^ dynamic  load  vector  that  involves  products  of  the  mass  matrix 
and  vectors  of  known  quantities  such  as  displacements,  velocities,  and  accelerations,  and 
At  and  P are  the  time  step  and  a parameter  arising  from  the  particular  Integrator  used, 
respectively. 

Unlike  an  explicit  integrator  there  is  a distinction  between  the  tangent  modulus  and 
initial  strain  approaches  when  using  on  implicit  Integrator.  This  distinction  is  in  the 
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same  sense  as  that  existing  for  static  problems.  That  is,  the  effect  of  plasticity  can 
enter  directly  into  the  stiffness  matrix  and  thereby  the  coefficient  matrix,  R,  or  it  can 
appear  as  an  effective  load  vector.  Equation  (IT)  represents  the  tangent  stiffneas  approach 
since  the  coefficient  matrix  explicitly  contains  the  tangent  stiffness  matrix  K.^,.  In  this 
case  the  matrix  R must  be  reassembled  and  a complete  solution  obtained  in  each  time  step. 

As  the  mesh  size  increases  this  becomes  increasingly  expensive  and  ultimately  dominates 
the  calculations. 

Using  the  initial  strain  approach  leads  to 

S:Au  - AP  + Q . ♦ AQ  + R (18) 

d 

'V  M 

where  K = K + ^ elastic  stiffness  matrix,  and  AQ  is  the  plastic  pseudoload 

vector.  If  a constant  time  step  approach  is  used  the  coefficient  matrix  remains  constant 
so  that  subsequent  solutions  require  only  forward  and  back  substitutions  of  triangular 
matrices.  This  procedure  is  competitive  with  an  explicit  technique  since  after  the 
first  step,  computations  involve  only  solutions  of  the  factored  coefficient  matrix,  the 
formation  of  the  dynamic  load  vector  and  the  formation  of  AQ.  It  must  be  noted  that 
the  time  step  used  in  this  procedure  must  be  controlled  by  the  linearization  implicit  in 
using  estimated  values  of  plastic  strain  increments,  rather  than  accuracy  consideration 
defined  entirely  by  the  integrator. 

An  inner  loop  iterative  procedure  can  be  developed  on  the  basis  of  Eqs.  (17)  or  (l6). 
Iterations  are  performed  until  the  structure  is  in  equilibrium  to  within  some  predetermined 
tolereuice.  This  tolerance  can  be  based  on  a measure  of  the  change  of  displacement  incre- 
ment [6Ti]  or  on  the  work  done  by  the  residual  load  vector  [63],  In  either  case,  if  the 
procedure  does  not  converge  within  a prescribed  number  of  Iterations  the  time  step  is  re- 
duced. Thus,  the  iterations  ceui  define  a variable  time  step  approach  In  which  the  criterion 
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for  reducing  the  time  step  Is  based  on  the  system  nonllnearltles.  It  Is  our  opinion  that 
accuracy  checks  should  also  be  included  based  on  the  current  system  dynamics  in  the  same 
manner  as  used  in  predictor-corrector  methods.  This  has  been  suggested  in  [591  but  to 
our  knowledge  has  not  been  implemented  and  tested  in  a general  purpose  program. 

There  are  currently  a number  of  implicit  operators  that  are  frequently  used  for 
structural  analysis.  Among  these  are  the  Newmark-6  family  (67]i  Wllson-0  [68],  Houbolt 
[69]  and  the  more  recently  developed  stiffly  stable  methods  by  Piurk  [70].  No  one  method 
at  this  point  seems  to  have  been  accepted  as  the  "best"  one  for  all  problems. 

SPECIAL  CONSIDERATIONS 

One  of  the  distinguishing  features  of  many  practical  problems  associated  with  plastic 
analyses  is  that  the  nonlinearities  ceui  be  contained  in  loceilized  regions.  Because  of 
this,  substructuring  in  a number  of  forms  cam  be  an  important  factor  in  reducing  com- 
putatloned.  costs.  The  most  straightforward  approach  is  to  simply  use  a Guyan  [71]  re- 


duction technique  to  reduce  the  number  of  degrees  of  freedom  in  the  entire  model,  partic- 
ularly in  the  elastic  region  and  then  to  recover  tresses  and  strains  only  in  the  limited 
region  where  plastic  flow  is  postulated.  Mixed  methods  [72,  73]  are  also  used,  where  modal 
techniques  are  utilized  in  some  well  defined  lastic  region  with  direct  integration  used  in 
the  plastic  region.  Other  methods  are  currently  being  developed  [7l*-77]  for  use  in  mixed 
media  problems  such  as  fluid  structure  interaction  and  hold  promise  for  dynamic  plastic 
emalysis.  Hughes  [7^]  has  demonstrated  a mixed  implicit-explicit  method  with  a predictor 
corrector  explicit  operator  and  a Newmark  family  implicit  integrator.  Within  our  context 
this  method  may  be  useful  for  elastic-plastic  dynamic  analysis  with  the  implicit  in- 
tegrator used  in  the  elastic  region  and  the  explicit  used  in  the  elastic-plastic  region. 
Other  techniques  that  allow  different  time  steps  in  different  regions  [75,  76,  77].  warrant 
investigation  since  simple  explicit  integration  with  time  steps  consistent  with  the  path 
dependence  attributes  of  plasticity  theory  could  be  combined  with  larger  time  step  implicit 
integration  in  a larger  elastic  region. 
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MODAL  METHODS 


Modal  superposition  for  linear  elastic  systems  nequlres  the  equations  of  motion  to  be 
transforswd  into  an  uncoupled  form  in  which  only  a portion  of  the  eigenvalues  and  eigen- 
vectors are  retained.  To  reduce  computational  costs  this  transformation  may  be  preceded 
by  a condensation  procedure  that  eliminates  a prescribed  set  of  degrees  of  freedom.  One 
of  two  procedures  is  usually  employed  to  eliminate  the  unwanted  information;  the  first  is 
the  Guyan  reduction  scheme  [7l]»  end  the  second  procedure  is  a static  (sero  mass)  conden- 
sation method.  Both  procedures  result  in  an  approximate  representation  of  the  mass  amd 
stiffness  characteristics  of  the  discrete  model.  As  in  the  case  of  substructuring  in  a 
static  analysis,  care  must  be  taken  to  avoid  eliminating  regions  where  plasticity  may 
develop.  The  effects  of  plasticity  may  be  treated  within  the  framework  of  the  tangent 
modulus  or  residual  force  method. 

The  use  of  modal  superposition  for  nonlinear  problems  appears,  at  first  glance,  to 
violate  the  well-known  fact  that  superposition  principles  are  not  applicable  to  nonlinear 
systems.  However,  nonlinear  dynamic  analysis  by  modal  superposition  requires  some  addi- 
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out  the  analysia.  Only  the  residual  force  due  to  plasticity  vould  have  to  be  transformed 
In  each  time  step. 

CLOSING  COMMENTS 

The  computatloned  capability  available  for  the  plastic  analysis  of  structures  has 
experienced  a tremendous  growth  during  the  past  decade.  Indeed,  the  level  of  structural 
ancdysls  capability  that  has  been  achieved  has  at  times  outstripped  our  ability  to  describe 
accurately  complex  material  behavior  such  as  cyclic,  time,  and  temperature  dependent  plas- 
ticity. Prior  to  the  development  of  the  programs  now  available,  the  designer  or  Malyst 
confronted  with  a problem  Involving  material  nonllnearltles  was  left  with  a choice  of 
using  his  engineering  Judgment  alone  or  in  conjunction  with  potentially  expensive  labora- 
tory tests.  He  now  has  the  f’orther  option  of  performing  numerical  emalysis  to  gain  insight 
into  the  behavior  of  the  structure. 

Most  nonlinear  analysis  programs,  with  the  exception  of  a few,  have  been  developed 
as  a spin-off  of  existing  programs  that  were  originally  designed  for  linear  structural 
analysis.  Although  this  development  is  a nature^  one,  this  added  dimension  of  generaLLlty 
has  placed  great  responsibilities  on  the  user  of  such  programs.  Perhaps  the  greatest  asset 
of  these  programs,  l.e.,  their  ability  to  solve  sophisticated  problems  also  represents  a 
potential  liability,  l.e.,  they  always  produce  numbers.  The  user  must  still  exercise  en- 
gineering Judgment  in  order  to  interpret  the  results  meaningfully.  Hopefully,  the  analytic 
results  will  confirm  these  feelings  and  provide  him  with  additional  insight.  However,  he 
now  has  the  luxury  of  having  his  intuition  fall  him  without  suffering  the  consequences  of 
a catastrophic  failure  or  an  overdesigned  system. 

For  the  analyst  and  researcher  of  fundamental  structural  and  material  phenomena,  the 
combination  of  advanced  numerical  analysis  procedures  and  the  appropriate  modeling  of  plas- 
ticity will  continue  to  be  used  to  gain  insight  into  the  significance  of  some  of  the  many 
factors  associated  with  plastic  deformation.  A perusal  through  the  preceedings  of  a sym- 
posium dedicated  to  the  numerical  modeling  of  manufactui'lng  processes  [T9]  and  [8o]  should 
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provide  the  reader  with  an  appreciation  of  the  application  of  some  of  the  models  and  pro- 
cedures discussed  here.  These  models  have  contributed  significantly  toward  the  realistic 
treatment  of  such  problems  as  welding  and  metal  forming. 
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Fig.  3 Piecewise  linear  hardening  (a)  Independent  loading  surfaces 
(b)  Independent  loading  surfaces  (c)  Special  case  of  (b) . 
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